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ABSTRACT 

In an effort to probe the globular cluster (GC) system of an isolated elliptical galaxy, 
a comprehensive analysis of the NGC 821 GC system was performed. New imaging 
from the WIYN Mini-Mosaic imager, supplemented with HST WFPC2 images reveals 
a GC system similar to those found in counterpart ellipticals located in high density 
environments. To put these results into the context of galaxy formation, a robustly- 
determined census of GC systems is presented and analysed for galaxies spanning a 
wide range of masses (>M*), morphologies and environments. 

Results from this meta-study: (1) confirm previous findings that the number of 
GCs normalized by host galaxy stellar mass increases with host stellar mass. Spiral 
galaxies in the sample show smaller relative GC numbers than those of massive ellip- 
ticals, suggesting the GC systems of massive ellipticals were not formed from major 
spiral-spiral mergers; (2) indicate that GC system numbers per unit galaxy baryon 
mass increases with host baryon mass and that GC formation efficiency may not 
be universal as previously thought; (3) suggest previously reported trends with en- 
vironment may be incorrect due to sample bias or the use of galaxy stellar masses 
to normalize GC numbers. Thus claims for environmentally dependent GC forma- 
tion efficiencies should be revisited; (4) in combination with weak lensing halo mass 
estimates, suggest that GCs formed in direct proportion to the halo mass; (5) are 
consistent with theoretical predictions whereby the local epoch of re-ionization did 
not vary significantly with environment or host galaxy type. 

Key words: galaxies: star clusters: globular clusters; individual: NGC 821; evolution 
& formation; elliptical and lenticular, cD; cosmology: early Universe 



1 INTRODUCTION 

Observational evidence suggests galaxies in high-density re- 
gions of the Universe began to form stars at earlier epochs 
(redshifts of 2: ~ 5 or ~ 11 Gyr) and finish within shorter 
time periods (of roughly a few Gyr) compared to galaxies 
of low-density environments. On average, galaxies in low- 
density environments started to form at 2: ~ 2 and sus- 

* Based upon data from the WIYN Observatory, which is a joint 
facility of the University of Wisconsin-Madison, Indiana Univer- 
sity, Yale University, and the National Optical Astronomy Ob- 
servatories. Also includes analysis of observations made with the 
Hubble Space Telescope obtained from the ESO/ST-ECF Science 
Archive Facility. 

f Email: lspitler@astro.swin.cdu.au 



tained low-level star formation rates even to the present day 
(Trager et al. 2000; Kuntschner et al. 2002; Proctor et al. 
2004; Thomas et al. 2005). What remains unclear are the ex- 
act conditions responsible for such an early divergence in the 
star formation histories of galaxies in high and low-density 
environments. Globular clusters (GCs) are among the oldest 
(> 10 Gyr) stellar systems in the Universe and their forma- 
tion is associated with major star formation epochs (see the 
recent extragalactic GC system review of Brodie & Strader 
2006). As such, systems of GCs provide unique observational 
tracers for investigating the environmental impact on galax- 
ies at the onset of their formation (West et al. 2004). 

Nearly all types of galaxies host GC systems that are bi- 
modal in their observed optical colour distributions (see e.g. 
Peng et al. 2006; Strader et al. 2006). Since GCs are mostly 
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old, these bimodal colour distributions are interpreted as 
evidence for bimodal metallicity distributions, which may 
imply two major star formation epochs took place in most 
galaxies before the Universe was ~ 3 Gyrs old. Constrain- 
ing the absolute age or the relative time between these two 
star formation epochs would provide invaluable constraints 
on galaxy formation, but is currently a formidable task (see 
e.g. Strader et al. 2005). 

An alternative method for probing the early epochs is 
to determine the relative numbers of GCs across different 
galaxies. This gives an indication of GC formation efficiency, 
which in turn helps characterise galaxy star formation histo- 
ries as a function of galaxy mass, environment, and morphol- 
ogy. The present work aims to characterise the GC system 
hosted by the nearby, isolated elliptical NGC 821. Obser- 
vational constraints gleaned from this GC system are an 
important contribution as nearly all previous work exam- 
ined GC systems in the high-density environments of galaxy 
groups and clusters. 

NGC 821 is an elliptical (E6) galaxy (Mv = -21.0; de 
Vaucouleurs et al. 1991) with a surface brightness fluctu- 
ation (SBF) distance modulus of 31.75 ± 0.17 (~ 22 Mpc; 
Tonry et al. 2001; corrected according to Jensen et al. 2003). 
It reveals a small edge-on disc both photometrically (Lauer 
1985) and kinematically (Emsellem et al. 2004) . Stellar ages 
range from 4.2 to 12.5 Gyrs in the literature, which Proctor 
et al. (2005) show to depend on the region of the galaxy 
sampled: the central regions are ~4 Gyrs old rising to 12.5 
Gyrs at 1 effective radius. Proctor et al. concluded that the 
central starburst event represents only a few percent of the 
total galaxy mass. Furthermore, they conflrm a steep metal- 
licity gradient and the alpha-element gradient was found to 
be flat with radius. No cold gas component has been detected 
in NGC 821 (Georgakakis et al. 2001; Nakanishi et al. 2007). 
Based on planetary nebulae kinematics, Romanowsky et al. 
(2003) claim NGC 821 has little, or no, dark matter out to 
~ 4 effective radii (c./. Dekel et al. 2005 and Douglas et al. 
2007). 

After a detailed photometric analysis of the NGC 821 
GC system, an in-depth comparison between this isolated 
GC system and those found in other environments is pre- 
sented. The implications for galaxy formation are discussed. 



2 OBSERVATIONS AND INITIAL DATA 
REDUCTIONS 

Imaging of NGC 821 was acquired over three nights in late 
2005 under good seeing conditions (average 0.7") with Uni- 
versity of Wisconsin-Madison time using the Mini-Mosaic 
(9.6' X 9.6') imager installed on the 3.5-meter Wisconsin 
Indiana Yale NOAO (WIYN) Telescope. The fihers used 
(B, V, and R) and exposure times of the final image prod- 
ucts are summarised in Table [T] For data reduction, the 
IRAF/PyRAFQ 

package MSCRED was used following the 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. PyRAF is a product of the 
Space Telescope Science Institute, which is operated by AURA for 
NASA. 



Table 1. Exposure times for a given camera, target and filter. 



AA/'TV NT 
VV i 1 IN 


WDS. uare. zuuo inov ou - uec z 


NGC 821 - B 
NGC 821 - V 
NGC 821 - R 


9 X 900s 

4 X 900s 

5 X 900s 


WFPC2 


HST Data: U307120*B 


NGC 821 - V 
NGC 821 - I 


4x350s + lxl20s 
6x230s + lx70s 


WFPC2 


HST Data: U4AN1[1-2]01B 


Background - V 
Background - I 


2 X 2400s 
2 X 2400s 




Figure 1. Digital Sky Survey image centred on NGC 821 show- 
ing the image footprints of the WIYN (square region) and two 
WFPC2 pointings (series of four boxes). The circle illustrates the 
spatial extend of the GC system estimated in i]3.5l 



WIYN Mini-Mosaic online reduction notes by S. Kafka. Fig- 
ure[T]shows the orientation and coverage of the single WIYN 
pointing. Observations were designed so that a 10th magni- 
tude (B-band) star fell in the 7.12" gap between the two 
Mini-Mosaic CCDs, thus minimising severe saturation ar- 
tifacts. The native pixel scale of 0.141" pixel"^ was main- 
tained for the final science images. 

An archive HST WFPC2 dataset in V and I bands was 
also analysed to complement the WIYN images where the 
light from NGC 821 and the bright star significantly de- 
creased object detection rates. As shown in Figure [T] this 
field is centred on NGC 821. A second WFPC2 pointing 
that only just overlaps with the WIYN image is shown 
in the same figure. This second dataset served as an esti- 
mate of the background contamination level for the central 
WFPC2 pointing. Although the second pointing exhibits a 
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fainter completeness limit than the WFPC2 pointing con- 
taining NGC 821, the only relevant analysis (characteris- 
ing the GC luminosity function, see H3.7P robustly accounts 
for this difference. The central WFPC2 image was down- 
loaded from the ESO HST Science Association data archive. 
The WFPC2 background pointing was produced manually 
with IRAF IMCOMBINE and designed to have a total ex- 
posure time similar to that of the central WFPC2 dataset. 
Table [T] provides the effective exposure times of the WFPC2 
datasets. 



3 DATA ANALYSIS 

3.1 Outline 

The primary objective of this work is to investigate GC for- 
mation efficiency in an isolated elliptical. More specifically, 
the number of NGC 821 GCs normalised by the mass of 
NGC 821 is a useful quantity to gauge GC formation efficien- 
cies between difl^erent galaxies. Accurately determining this 
value requires a robust estimate of the NGC 821 GC pop- 
ulation, motivating the analytical procedure outlined here 
and carried out in the remainder of 331 

Most GCs at the distance of NGC 821 will appear as 
point sources on everything but the small region covered by 
the high-resolution WFPC2-PC images and will therefore 
resemble Galactic stars and unresolved galaxies. An objec- 
tive selection criteria is therefore needed to remove these 
contaminates from the raw object dataset. This proceedure 
will be described in tj3.2l and H3.4I Estimating the amount 
of contamination inevitably missed in this selection process 
is an important component to the GC number estimate (see 

MM- 

Observational biases need to be quantified and consid- 
ered during analysis. For instance, the detection rates vary 
significantly with galactocentric distance due to the lumi- 
nous centre of NGC 821 f H3.3p . Accounting for this motivates 
radial-dependant corrections to the GC luminosity function 
(GCLF; see tj3.7|l and the use of the central IfST imaging. 
To estimate the total number of metal-poor and metal-rich 
GCs, the varying proportion of these GC types again with 
galactocentric distance must be understood as well as any 
selection biases on colour ( tj3.6|l . 

This analysis comes together in i]3.8l and where robust 
population estimates of the NGC 821 GC system are pre- 
sented. Implications from this work are explored in 311 

3.2 GC Photometry 

3.2.1 WIYN 

Sources on the WIYN images were detected on each of the 
three bands after a median filter (of 31 x 31 pixel window) 
image was subtracted from the original images, thereby re- 
moving the diffuse light of NGC 821. Objects 3a above the 
global background were catalogued with the find task in 
DAOPHOTII (Stetson 1992). Flux from each object was ex- 
tracted with DAOPHOTII phot in an aperture designed to 
keep the subsequent aperture corrections small. Images still 
containing the diffuse NGC 821 light were used for the pho- 
tometry. Local background flux levels for individual objects 
was subtracted from the extracted flux and estimated from 



Table 2. WIYN photometric corrections for the finite apertures. 
Galactic extinction and atmospheric reddening, each in units of 
magnitudes. For each filter, the average airmass during observa- 
tions was used for such corrections. 



Filter 


Aperture 


Extinction 


Airmass 


B 


-0.14 ±0.04 


-0.48 


-0.22 


V 


-0.13 ±0.04 


-0.37 


-0.14 


R 


-0.19 ±0.03 


-0.30 


-0.08 



a 10-pixel wide annulus with an inner radius of 15 pixels (or 
~ 2.1") from the object centre. Aperture corrections were 
derived from point sources on the images. Corrections for 
Galactic foreground (with a reddening excess of E(B — V) 
= 0.11; Schlegel et al. 1998) and airmass extinction were 
also incorporated into the total flux measurements. Photo- 
metric zeropoints derived from standard Landolt (1992) field 
star observations were applied to the fiuxes yielding the final 
magnitudes values on the standard Johnson- Cousin photo- 
metric system. Photometric corrections are summarised in 
Tabled 

3.2.2 WFPC2 

Both WFPC2 datasets (see Table[l| were analysed following 
the aperture photometry procedure described in Larsen et 
al. (2001). To this end, a generalised WFPC2 data reduction 
pipeline provided by S. Larsen was employed. This pipeline 
makes use of the IRAF APPHOT package. Objects devi- 
ating 3cr from the background noise level were found and 
photometry was measured in the central WFPC2 field from 
a median filter (15 x 15 pixel window) subtracted image that 
effectively removed the diffuse light gradient from NGC 821. 
For the background pointing, individual background values 
were computed for each object in annuli (aperture sizes as 
in Larsen et al. 2001) centred on each object. All fluxes 
were calibrated to the standard photometric system accord- 
ing to the Holtzman et al. (1995) procedure and corrected 
for charge-transfer efficiency as presented in Stetson (1998). 

Corrections for Galactic foreground extinction were ap- 
plied with the same reddening excess value used for the 
WIYN photometry. Aperture corrections to the 0".5 Holtz- 
man et al. (1995) reference aperture were directly measured 
from the objects passing GC selection criteria (see 8 ^3.41 for 
this criteria). The median magnitude difference between the 
extraction and reference aperture of these candidate GCs 
on the WFPC2-WF CCDs was added to all fluxes to cor- 
rect for light found beyond the extraction radius. Separate 
aperture corrections for the WFPC2-PC CCDs were neces- 
sary because of its different point spread function. There are 
not enough GC-like objects in the background WFPC2-PC 
CCD, thus a full extraction aperture of 0".5 was used for 
the 2 objects found within the WFPC2-PC field. WFPC2 
photometric corrections are summarised in Table O 

3.3 Completeness Tests 

Completeness tests were performed to quantify the magni- 
tude limits of each dataset. Using the IRAF DAOPHOT 
package, artificial objects produced from empirical point 
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Table 3. WFPC2 photometric corrections for finite apertures 
and Galactic extinction. 



Target 


WFPC2 


Filter 


Aperture 


Extinct. 




Chip 




Correction 




NGC 821 


PC 


V 


-0.10 


-0.36 


NGC 821 


WF 


V 


-0.01 


-0.36 


NGC 821 


PC 


I 


-0.17 


-0.22 


NGC 821 


WF 


I 


-0.19 


-0.22 


Background 


WF 


V 


-0.23 


-0.36 


Background 


WF 


I 


-0.25 


-0.22 



spread functions were placed on to the original science im- 
ages in such a way that crowding from existing and other 
artifical objects could not influence the test results. For each 
0.2 magnitude step between the expected magnitude range 
of GCs, 15000 and 500 objects of the same magnitude were 
added to the WIYN images and each WFPC2 CCD, re- 
spectively. Since the WFPC2 data pipeline requires both a 
V and I-band image as input, the colour of the artificial 
objects was taken to be V— I = 1.1. This colour is typi- 
cal of GCs (e.g. Kundu & Whitmore 2001 and Larsen et 
al. 2001) The appropriate photometry extraction procedure 
(for WIYN see ^TTTl WFPC2 Sl3X2j were carried out on 
each artificial field. Completeness tests are given in Figure (2] 
When appropriate, variations with galactocentric radius in 
the completeness function are considered in the subsequent 
analysis. 

3.4 Globular Cluster Selection Criteria 

In the following two subsections, criteria are defined that 
effectively reduce the WIYN and WFPC2 raw object list 
to mostly GC candidates. Part of the final GC catalogue is 
given in Tabled the full version of which is available from 
the authors or available online at the CDS. For objects in 
common between the WIYN and WFPC2 datasets, the tri- 
band WIYN selection took precedence over the WFPC2 GC 
V— I selection, unless the object was 0.'4 from NGC 821 in 
which case the WFPC2 V— I selection was used. Aside from 
the data presented in Table [3] and the analysis where merg- 
ing was necessary (e.g. determining GC subpopulation pro- 
portions, tj3.6l and the spatial distributions, i)4.ip . for the 
subsequent NGC 821 GC analysis, the WIYN and WFPC2 
datasets were analysed independently of each other. 

3.4.1 WIYN 

Spurious detections such as image artifacts and statistical 
events were removed when individual objects were matched 
by their positions between the B, V and R-band WIYN 
images. An extended-source culling of the remaining 1073 
objects was carried out by removing objects that show an 
excess of light beyond the extraction aperture compared to 
point-sources (e.g. Holtzman et al. 1996). This involves an 
analysis of the magnitude difference between the adopted 
extraction aperture and the outer aperture that defined the 
aperture corrections. Figure [3] shows sequences of relatively 
bright point sources (most likely stars) that make up the 
fiducial point-source magnitude difference and define the 




20 22 24 26 

V 



Figure 2. Completeness functions from tests described in the 
text. Top panel: Results for the WIYN B-, V- and R-band images, 
represented by solid, dotted and dashed lines, respectively. To 
facilitate this comparison, the B and R-band magnitudes were 
converted to V magnitudes with the mean GC colours of B— V 
= 0.69 and V— R = 0.45. Object recovery on the B-band image 
approaches the 50% level for objects ~ 1 mag. brighter than the 
V and R-bands. Bottom panel: WFPC2 completeness test results 
for the central and background pointings are represented by solid 
and dotted lines, respectively. 



aperture corrections. Measurement uncertainties were ac- 
counted for by using a magnitude difference selection that 
depends on the final magnitude of the object. The FWHM 
difference between point and extended sources suggests the 
size selection employed will remove any object at the dis- 
tance of NGC 821 with a effective radius greater than ~ 9 pc. 
Less than 4% of old Milky Way show sizes greater than 9 pc. 
A total of 746 point-sources remain after the application of 
this aperture difference selection criteria. 

A two-colour GC selection technique was used to further 
constrain the GC candidate list (e.g. Rhode & Zepf 2001) 
to the same region occupied by known GCs in B— V versus 
V— R colour space. Rhode & Zepf (2001) used Galactic GCs 
as their basis sample. Milky Way GCs have a relative defi- 
ciency of metal-rich GCs, thus the WIYN sample is instead 
compared to accurate HST photometry of 659 NGC 4594 
GCs (Spitler et al. 2006), which are partly resolved in the 
HST images. The NGC 821 and NGC 4594 datasets are 
presented in Figure |3] A linear relation was defined from 
the NGC 4594 GC colour-colour information and the final 
NGC 821 GC candidates are those with V— R colours falling 
within 2a (where a is the intrinsic scatter of NGC 4594 GCs) 
of this relation. B— V colour limits of 0.55 and 1.05 were also 
imposed and correspond to [Fe/H] values of —2.5 and -fO.2 
for Galactic GCs, according to a linear fit to the data in the 
Harris (1996) catalogue. Following Rhode & Zepf (2001), if 
an object's colour deviated from the criteria by an amount 
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Table 4. Sample of the catalogue of astrometry and photometry for GC candidates detected around NGC 821 on the WIYN and central 
WFPC2 images. The full table is available at the CDS. Note those objects beyond a galactocentric radius of 4' from NGC 821 made up 
a contamination sample for the WIYN analysis (see i|3.5l l. 



ID 


RA (2000.0) 


Dec (2000.0) 


Rgc(arcmin) 


B 


crB 


V 


crV 


R 


aR 


/ 


al 


1 


02:08:21.425 


11:05:14.18 


5.53 


21.163 


0.004 


20.618 


0.004 


20.232 


0.003 






2 


02:08:20.657 


10:58:55.60 


0.79 


21.741 


0.012 


20.764 


0.006 


20.312 


0.010 


19.618 


0.005 


3 


02:08:25.373 


10:59:53.71 


1.06 


21.806 


0.007 


21.096 


0.007 


20.633 


0.006 








0.8 



0.6 



a 

I 

> 



0.4 



0.2 




0.5 1 
B-V 



1.5 



Figure 3. Selection criteria applied to remove extended sources 
(e.g. background galaxies) detected on the NGC 821 WIYN im- 
ages. The final B {top panel), V (middle), and R-band magnitudes 
{bottom) are plotted against the respective magnitude difference 
between the adopted extraction aperture and the one used to 
estimate the aperture corrections. In each panel, those objects 
with the largest differences are considered extended sources and 
are represented as open circles, while point-source candidates are 
shown as solid circles. In the online edition, the blue, green and 
red coloured open circles represent the flagged objects from ex- 
tended object selection in the B, V and R-band photometry, re- 
spectively. 



Figure 4. NGC 821 WIYN object V-R versus B-V plot. Solid 
circles represent the final GC candidates, which fall within 2a 
of the colour region defined by NGC 4594 GC HST photometry 
(small grey circles, Spitler et al. 2006). This criteria was relaxed 
by an amount equal to an object's colour error, as described in 
the text. Open triangles are WIYN objects failing these criteria 
and are probable stars or galaxies. 



Table 5. Ultra-compact dwarf galaxy candidates. The fourth 
entry was found on a WFPC2-WF chip, thus there is a 
chance the object is partially resolved. Size measurements with 
IRAF/IMEXAM indicates no difference between the candidate 
and two bright stars on the same chip. This puts a ~ 5pc upper- 
limit on the effective radius of this UCD candidate. 



equal to its formal colour error, it was still considered a GC 
candidate. 

Finally, objects were considered GC candidates if their 
V-band magnitude is fainter than V = 20 mags (i.e. Mv = 
— 11.75). This magnitude limit approximately corresponds 
to the boundary between GCs and more massive compact 
stellar objects (e.g. Hasegan et al. 2005; Mieske et al. 2006). 
It is noted that three (plus one on the WFPC2 central point- 
ing) point-sources with GC-like colours are brighter than V 
= 20 mags. These are candidate ultra-compact dwarf galax- 
ies and are presented in Table [S] Objects with photometric 
errors greater than 0.15 mags in any one band were also dis- 
carded, leaving 237 GC candidates from the WIYN images. 





R.A. 


Dec. 


V 


B-V V-I 


02 


08:05.135 


-1-10:57:35.13 


19.47 


0.96 


02 


08:34.641 


-1-11:00:38.03 


18.94 


0.54 


02 


08:16.697 


-1-11:01:24.51 


17.93 


1.07 


02 


08:17.685 


-1-11:00:31.84 


19.60 


1.20 



3.4.2 WFPC2 

The adopted GC selection criteria for the two WFPC2 fields 
are identical. Object with elongated or extended structure 
(determined by eye and with IRAF/IMEXAM) were culled. 
Remaining objects detected on the WFPC2 fields are shown 
in the colour- magnitude diagram in Figure [5] GC colour se- 
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Figure 5. Colour-magnitude diagram of objects detected in the 
HST WFPC2 pointings. Solid circles and open triangles repre- 
sent point sources found in the central NGC 821 and background 
pointings, respectively. GC selection limits are indicated with ver- 
tical lines at V— I = 0.7 and = 1.3 and a horizontal line at V 
= 20.0. 



Figure 6. Raw WIYN NGC 821 GC surface density profiles. The 
bottom panel shows a linear representation of the data and the 
top panel is in logarithmic form. The consistency between the 
outer 6 points suggests the data reaches a constant background 
level by Rgc ~ 4 arcmin (= 26 kpc; indicated by the dashed 
vertical line). Uncertainties are from Poisson statistics. 



lection of 0.7 < V— I ^ 1.3 mags were applied and corre- 
sponds to roughly the same metallicity range as the WIYN 
B~V colour limits. One object meeting these requirements 
was brighter than V = 20 mags and was removed from fur- 
ther analysis, but is recognised as a candidate ultra-compact 
dwarf galaxy (see Table [SJ. A total of 156 and 65 GC candi- 
dates remain after these criteria were applied in the central 
NGC 821 and background WFPC2 pointings, respectively. 



the GC selection criteria in this pointing conveniently define 
the contamination level for the central NGC 821 WFPC2 
pointing. A combination of differing exposure times in these 
two datasets and the presence of NGC 821 in the central 
pointing undoubtedly leads to differing detection rates. This 
is only relevant to characterising the GCLF (see i]3.7|l and 
so the GCLF analysis employed robustly accounts for this 
issue. 



3.5 Contamination Estimates 

The amount of contamination in the WIYN NGC 821 
dataset was estimated from the objects found in the region 
beyond the spatial extent of the GC system. Figure [6] shows 
the surface density profile from the raw counts of GC can- 
didates in 12 radial annuli of equal widths (0'.5) centred on 
NGC 821. Image boundaries in each band and masked re- 
gions where objects could not be detected (e.g. bright stars) 
were excluded in the annulus areal coverage calculations. It 
is apparent in the figure that a constant background level of 
1.5 — 2.0 GC candidates per square arcmin begins at a pro- 
jected distance of ~ 4 arcmin (4' = 26 kpc) from NGC 821. 
This galactocentric radius (Rgc) is taken as the "edge" of 
the GC system, despite the probable existence of a few GCs 
beyond this nominal edge. Figure [1] demonstrates that the 
GC system is effectively covered by the WIYN image, thus 
no further areal corrections are necessary for the GC popula- 
tion estimates. The adopted GC system and contamination 
unmasked regions have total areal coverage of 48.6 and 41.4 
square arcmin, respectively. 

As shown in Figure[l] the background WFPC2 pointing 
is also beyond the GC system edge, thus all objects passing 



3.6 Colour Distributions 

Figure [7] shows colour histograms of NGC 821 GC candi- 
dates in the WIYN and WFPC2 pointings. GC colour bi- 
modality is most apparent in the central WFPC2 pointing, 
while the WIYN colour distribution shows only a slight in- 
crease over the estimated contamination level at colours con- 
sistent with metal-rich GCs. 

A KMM test for bimodality (Ashman et al. 1994) sug- 
gests a marginal probability of bimodality in the WFPC2 
distribution (with a p-statistic of 0.0651) and no statisti- 
cally significant evidence (p-statistic of 0.6672) for bimodal- 
ity in the WIYN B— R distribution. Two important facts 
must be considered before interpreting these KMM results; 
the established observational property of GC systems that 
red (metal-rich) GCs tend to be more centrally concen- 
trated (e.g. Geisler et al. 1996) and a poor detection rate 
near NGC 821 in the WIYN images. The first point leads 
to a detection bias towards red GCs on the WFPC2 im- 
age, while the combination of the two biases the WIYN de- 
tections against red GCs. Overcoming these biases requires 
combining the WIYN and WFPC2 GC colour information, 
thereby using the WFPC2 data to "fill-in" the central region 
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Figure 7. GC and background contamination colour histograms. 
Top panel: Shows the B— R histograms of candidate GCs (solid 
lines) and contamination objects (dotted lines) extracted from the 
WIYN pointing. Bottom panel: V— I distribution from the central 
WFPC2 pointing and the contamination level estimated from the 
WFPC2 background image. The split between the blue and red 
GCs occurs at B— R ~ 1.3 and V— I ~ 1.1. Blue GCs dominate 
the WIYN distribution, while the proportion of red GCs is nearly 
half the total in the central NGC 821 WFPC2 dataset. 



rejected at a confidence level of 99.998% (with a p-statistic of 
0.00160), strongly supporting the apparent V— I bimodality 
in the figure. The KMM estimate for the two colour peaks is 
V— I = 0.99 and = 1.18, where the colour peak uncertainties 
are likely dominated by the zeropoint errors associated with 
the conversion from B-R to V-I colours: 0.14 mag. 

According to the KMM test, approximately 30% of the 
total falls into the red GC subpopulation. Figure |8] also 
shows the combined background contamination datasets, 
which similarly restricted in magnitude-space to account for 
the two biases detailed above. The red (V— I > 1.1) contam- 
ination level corresponds to approximately half the number 
of red GCs in the Figure, while the fraction of estimated 
contamination among the blue GCs is significantly lower. 
This perhaps indicates the proportion of red GCs is some- 
what smaller than the 30% derived with the KMM test, 
suggesting this value is actually an upper limit. 

To test the robustness of these KMM results as they 
relate to the uncertainty associated with the B— R to V— I 
colour conversion (see Eq. 1), a 500-run Monte Carlo sim- 
ulation was performed. The WIYN objects were converted 
to V—I colours using 500 different linear relations derived 
from randomly deviating the zeropoint and slope of Equa- 
tion 1 by their formal errors. According to a KMM analysis 
of the simulated datasets, a strong rejection of unimodality 
bimodality (i.e. a KMM p-statistic < 0.05; Ashman et al. 
1994) is found in ~ 80% of the runs, marginal evidence for 
bimodality occurs in 5% and the remaining 15% show no 
bimodality. Thus, a strong case for colour bimodality in the 
NGC 821 GC system is supported at a la confidence level. 



where detections fall on the WIYN image, near NGC 821. 
The WIYN B— R GC colours were therefore converted to 
the WFPC2 V—I colours with an empirical conversion from 
GC candidates in common between the datasets: 

V-7 = 0. 24 ±0.14 + 0. 65 ±0.11(B-i?) (1) 

Objects in common between the WIYN and WFPC2 
datasets were only counted once. 

Two additional observational biases must also be ac- 
counted for when robustly characterising the GC colour dis- 
tribution. First, the combined colour data were corrected for 
colour biases arising because a single band limits the sam- 
ple, which spans a range of colours. The following faint lim- 
its therefore need to be enforced so that GCs of all colours 
were detected at the > 95% completeness level: WIYN R 
< 23.0 mag and WFPC2 V < 23.5. Such a restriction also 
reduces uncertainties associated with contamination, which 
is generally faint. Second, metal-rich GCs have higher mass- 
to-light ratios due to emission line blanketing, so these GCs 
should be fainter (by B — 0.3 and R = 0.2 mags; Ashman 
et al. 1995) compared to a metal-poor GC of identical mass. 
Colour-dependent magnitude selections were therefore ap- 
plied: VvKFPC2 < 3/2 V-I +21.8 and Rwiyn < 2/3 B-R 
+22.4. 

The combined GC data, corrected for these biases (to- 
talling 87 objects), is presented in Figure [S] A bimodal 
colour distribution is apparent: a large, "blue", peak in the 
histogram is clearly observed at V—I ~ 1 and a smaller "red" 
peak is centred at V—I ~ 1.15. Results from KMM analysis 
on this combined colour dataset suggest unimodality can be 



3.7 Globular Cluster Luminosity Function 

The mean magnitude of a particular GC system is fairly 
constant for all galaxies for which it has been observed (e.g. 
Harris 2001; Richtler 2003; Jordan et al. 2007b). The shape 
of the GCLF can be approximated with a Gaussian distri- 
bution and in principle, the Gaussian "peak" or turnover 
magnitude (TOM) of a GCLF provides a useful quantity 
for comparison between GC systems. Richtler (2003) aggre- 
gated WFPC2 GCLF TOM values and produced a "univer- 
sal" TOM of Mv = —7.35 (here adjusted for the Jensen et 
al. 2003 SBF correction) with an observed RMS of 0.25 mag. 
Observational evidence suggests the GCLF does not vary 
with Rgc (see refs. in de Grijis & Parmentier 2007), thus 
the GCLF in the central WFPC2 region should be similar 
in form to the GCLF derived from the surrounding WIYN 
coverage. 

Seeker & Harris (1993) designed a maximum-likelihood 
technique to characterise the GCLF. This program corrects 
for magnitude incompleteness, contamination, and compen- 
sates for a different completeness function for the contami- 
nation dataset. It produces a best-fitting GCLF TOM and 
dispersion assuming either a Gaussian or Student ts (is) dis- 
tribution. Completeness functions and contamination levels 
differ between the WIYN and WFPC2 datasets, thus the 
analysis in this section makes use of the separate GC candi- 
date lists from the respective images. 

Assuming the Richtler (2003) universal TOM applies to 
NGC 821, GCs at the TOM magnitude should show an ap- 
parent magnitude of V~ 24.4 + 0.3. Figure [9] shows that the 
completeness level at this magnitude is 45% and 75% for the 
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Figure 8. Histogram of combined WIYN and WFPC2 V-I 
colour values. Solid and dotted histograms are for the GC can- 
didates and objects in the background contamination fields, re- 
spectively. See text for magnitude limits imposed to remove colour 
biases. WIYN B— R photometry was converted to WFPC2 V— I 
with an empirical relation derived from common objects. Output 
from a KMM analysis of this data strongly supports the apparent 
bimodality and indicates ~ 30% of the GCs in NGC 821 can be 
associated with the red subpopulation. A Monte Carlo simulation 
supports GC colour bimodality at the lo" level. 



WIYN, and WFPC2 datasets, respectively. Because only a 
fraction of the NGC 821 GCs at the expected TOM are ac- 
tually observed, the uncertainty associated with the GCLF 
parameters from the Seeker & Harris program increases by 
a non-negligible amount, as described in Seeker & Harris 
(1993). Nevertheless, applying this program to the WFPC2 
GC candidate and contamination datasets, assuming a 
form of the GCLF (which provides a better fit than a Gaus- 
sian), yields the following parameters: Vtom = 24.2 ± 0.3 
and a dispersion of 0-45 — 1.0 ± 0.2. A Gaussian fit produced 
the following values: Vtom = 24.3 ± 0.4 and ct = 1.2 ± 0.2. 
The best-fit TOM values are consistent with the universal 
TOM. 

The Seeker & Harris (1993) program also estimates the 
number of GCs in the corrected, best-fit distribution. The 
WFPC2 ts fit corresponds to Ngc = 130 ± 20 withm the 
WFPC2 field of view. Estimating the total number of GCs 
in NGC 821 requires the extension of this analysis to larger 
galactocentric radii and hence the use of WIYN GC candi- 
dates. Unfortunately, because the WIYN 50% completeness 
limit approximately corresponds to the estimated GCLF 
TOM, the solution of the Seeker & Harris algorithm is much 
too uncertain for any practical use. Thus, a modified ver- 
sion of the Seeker & Harris (1993) program was designed 
to produce a GC number estimate for a specified GCLF 
TOM and dispersion. GCLF parameters derived from the 
WFPC2 data were assumed and applied to the WIYN data 
as described in the following paragraph. 



Figure 9. V-band NGC 821 GC candidate magnitude distribu- 
tion. For each panel, the left axis is for the histograms and the 
right is for the completeness functions (dashed lines). WIYN and 
WFPC2 datasets are presented in the top and bottom panels, 
respectively. Solid histograms are GC candidates and dotted his- 
tograms are the contamination datasets. The best-fit (5 distribu- 
tion is shown over the WFPC2 histogram. 



The B-band image limits the WIYN photometry (see 
i]3.3|l . thus the input completeness function to the Seeker 
& Harris program was the B-band completeness function 
converted to V with the mean GC colour: B— V — 0.69. 
Since the number of GCs and the WIYN completeness func- 
tion both vary significantly with Rgc, the WIYN GCLF fits 
were performed on the data in 0'. 5- wide annuli. Appropri- 
ate completeness functions are used for each Rgc interval and 
the final GC surface density estimates are given in Table |6l 
Similar analysis of the WFPC2 data yielded the two inner- 
most surface density values. Surface densities errors are from 
Poisson statistics added in quadrature with the uncertainty 
in GC number estimate. The latter value is dominated by 
the measurement uncertainty of the TOM, which is easily 
translated to a GC number error with the modified Seeker 
& Harris program. 



3.8 GC Radial Profile and Population Estimates 

The corrected NGC 821 GC surface density profile from data 
values in Table |S] is presented in Figure 1101 To derive an 
estimate of the GC population total that accounts for com- 
pleteness fraction variations with galactocentric radius, a 
function is fit to the surface density data, which is then inte- 
grated with respect to the area covered by the NGC 821 GC 
system. Three functions are fitted: a power-law {ace ocR"), 
de Vaucouleurs profile (ctgc = 10""*=^^^ '^+I3d^v^ ^nd King 
(1962) profile {aac = afe,ng[l/(l + r^/Plr^aV^^ - 1/(1 + 
■|^g2-ji/2j2-j ^ glance at the fits, presented in Figure [Till gives 
the impression that the King profile best represents the GC 
number density profile. However, both the de Vaucouleur 
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Table 6. Robust GC number density estimates from GOLF fits 
and the relevant masked image areas. Radii are the projected 
mean galactocentric distance from NGC 821 for a given annulus. 



Radius 
(arcmin) 



(arcmin"^) 



Unmasked Area 
(arcmin^) 



0.2 
1.1 
1.3 
1.8 
2.3 
2.8 
3.3 
3.8 



144. ± 38. 
20.4 dz 10. 
11.9 ±3.1 
6.40 ± 1.8 
6.78 ± 1.8 
2.00 ±0.8 
0.50 ±0.5 
0.00 ±0.4 



0.32 
4.74 
3.45 
5.31 
6.94 
8.51 
10.1 
11.7 
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Figure 10. GC surface density corrected for contamination 
and incompleteness. Uncorrected data points are given as open 
squares. Circles are the corrected surface density from WIYN 
data, except the two inner-most solid points are from the WFPC2 
data. The open circle is not fit due to the high incompleteness 
fraction near NGC 821 on the WIYN image. The outer six points 
define the contamination level. Formal errors are from Poisson 
statistics and GCLF measurement uncertainties added in quadra- 
ture. Best-fit power-law, de Vaucouleurs and King (1962) func- 
tions are given as solid, dotted and dashed lines, respectively. 



and King fits yield consistent GC population estimates when 
surface integrated to tlie GC system's edge: 330 and 335, re- 
spectively. The former profile is used in the subsequent anal- 
ysis since results derived from it are directly comparable to 
past studies (e.g. Rhode & Zepf 2004). 

To test the significance of the de Vaucouleur number es- 
timate, a Monte Carlo simulation was performed where the 
observed surface density points were randomly deviated by 
their formal errors (assuming Gaussian errors). The result- 
ing values were fitted, and the new profile was integrated 
as before. Repeating this 5000 times resulted in a distribu- 
tion of population numbers. The mean of the distribution 
is taken to be the final GC population total: N= 320 ± 45, 



where the quoted uncertainty is the standard deviation of 
the simulated number distribution. The present NGC 821 
GC system total is smaller than the Kundu & Whitmore 
(2000) number (N= 395 ± 94) they estimated within the 
central WFPC2 field of view. 

Analysis of the GC colour distribution (JSHJ suggests 
30% of the GCs are red, thus the NGC 821 GC system con- 
tains approximately 95 ± 15 red and 225 ± 30 blue GCs. 



4 RESULTS AND DISCUSSION 

4.1 The GC System of NGC 821 

GC specific frequency (Sjv) has been used for comparing 
the abundance of GCs across different galaxies (Harris & van 
den Bergh 1981) and is defined in Table[7] The NGC 821 GC 
system population estimated in i]3.8l (Nr;f7 — 320±45) corre- 
sponds to a V-band specific frequency of S]v = 1.32 ± 0.15. 
This value is slightly lower compared to values tradition- 
ally associated with low-luminosity elliptical galaxies (e.g. 
Harris 2001). However, recent efforts have yielded generally 
smaller S n values when a careful accounting for contamina- 
tion is performed (e.g. Rhode & Zepf 2004). Indeed, it will 
be shown in tj4.3l that the NGC 821 stellar mass-normalised 
number of GC associated is roughly consistent with the same 
values for similar-mass elliptical and spiral galaxies. 

The NGC 821 GC colour distribution shows a hint of the 
classical bimodal form, with a fairly typical fraction (30%) 
of red GCs. The subpopulation mean colour values (V— I 
— 0.97 and 1.18 mags.; errors ~ 0.14) are roughly consis- 
tent with the same values in galaxies of similar luminosities 
(Larsen et al. 2001; Strader, Brodie & Forbes 2004). While 
a young (~ 1 — 4 Gyr) central stellar population in NGC 821 
was confirmed by Proctor et al. (2005), they demonstrated 
this component only refiects ^ 10 per cent of the mass in 
the central regions of NGC 821. Beyond one effective ra- 
dius. Proctor et al. (2005) found consistently old ages (~ 12 
Gyr) from the integrated spectrum of NGC 821. With no 
evidence for a major star formation event in 12 Gyr, it is 
likely the majority of NGC 821 GCs are at least as old, as is 
believed to be the case for most GC systems (see references 
in Brodie & Strader 2006). The GC colour distribution can 
therefore directly be interpreted as evidence for a bimodal 
metallicity distribution, here shown with marginal certainty 
in NGC 821 for the first time (see Fig. O. 

Red GCs in Figure[TT]tend to be more centrally concen- 
trated, as is commonly found in GC systems. The GC sys- 
tem extends to a projected galactocentric radius of R^c ~ 4 
arcmin or ~ 26 kpc, which is a factor of ~ 5 greater than 
the effective radius of NGC 821 (50"; de Vaucouleurs et al. 
1991). 

In the region north-west of NGC 821 (see Figure 1111 
along the galaxy minor-axis) there is an apparent under- 
density of GC candidates relative to other azimuthal direc- 
tions. While the 10th magnitude foreground star (also just 
north-west of NGC 821) significantly influences detection 
rates on the WIYN images, the near-normal detection rates 
on the WFPC2 images sufficiently compensated for this ru- 
ined WIYN region. Thus, with no signiflcant azimuthal bias, 
the low numbers of GCs perpendicular to the major-axis of 
NGC 821 must be treated as a real trend, although small 
number effects are not yet ruled out. 
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Figure 11. WIYN colour composite image showing the locations of GC candidates (blue and red circles). The central WFPC2 image 
coverage is given. The image is oriented so that the major-axis of NGC 821 is parallel with the width of the page. A deficiency of GCs is 
apparent along the galaxy's minor-axis, which does not result from selection effects due to the bright foreground star immediately below 
the figure centre. 



Folded position angles (PAs) of individual GCs in the 
whole system and two subpopulations are presented in Fig- 
ure[l2l A PA of 0° and 90° (and -90°) corresponds to direc- 
tions along the galaxy major and minor-axis, respectively. 
Data in Figure [T^] confirms the spatial distribution of GCs 
in Figure [11] and apparent like orientation between the GC 
system and major-axis of NGC 821. This is consistent with 
previous observations (see summary in Brodie & Strader 
2006) where GC systems generally follow the shape of the 
spheroidal component of their host galaxy. For NGC 821 



the situation is likely more complex, as the dearth of GCs is 
most apparent in the direction north-west of NGC 821, and 
not so obvious to the south-east. It should be noted, there 
is also a hint of a major-axis orientation among the plane- 
tary nebula analysed by Romanowsky et al. (2003; see their 
fig. S2) . GC kinematic information would provide a valuable 
extra dimension to this analysis. 

Information gleaned from the NGC 821 GC system is 
generally consistent with a scenario where this galaxy pas- 
sively evolved since the formation of the bulk of the field 
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Figure 12. Folded position angle histograms of individual GCs 
associated with the whole system {top panel), blue {middle), and 
red {bottom) subpopulations. Errors are Poisson statistics. A PA 
of 0° corresponds to the direction along the galaxy's major-axis. 
A deficiency of GCs along the minor-axis is observed. 



stars (Proctor et al. 2005) and GCs roughly 12 Gyr ago. 
More specifically: a relatively steep radial profile is inconsis- 
tent with a dissipation-less major merger (Bekki & Forbes 
2006) and evidence for a standard bimodal GC colour dis- 
tribution is consistent with no late (< 5 Gyr) major GC for- 
mation epochs. Mild nuclear activity (Pellegrini et al. 2007) , 
negligible cold gas components (from limits on HI content 
and no CO emission detections: Georgakakis et al. 2001; 
Nakanishi et al. 2007) further support the idea that only 
small amounts of star formation took place during the re- 
cent past of NGC 821 (Proctor et al. 2005). An extended, 
but relatively subdued formation history is consistent with 
the conclusions from a spatially-resolved spectroscopic anal- 
ysis of 12 isolated early-type galaxies by Reda et al. (2007) . 



4.2 A New Compilation of Globular Cluster 
Systems 

To probe the relationship between GC system formation and 
environment a literature search for reliable GC system prop- 
erties was conducted. The compilations of Rhode, Zepf & 
Santos (2005) (hereafter referred to as RZS05) and Brodie 
& Strader (2006) served as the primary, but not exclusive 
source of literature references. Inclusion criteria are the same 
as in RZS05. Included GC system studies have: data for 50% 
of the GC system radial extent, observed in at least two fil- 
ters, total GC number and subpopulation proportion given 
or derivable, and a GC number estimate error less than 40%. 
Exceptions to the radial coverage restriction were given to 
NGC 5128 and M87, as these pose difficult observational 
targets, but were analysed sufficiently well (see Harris et 
al. 2006 and Tamura et al. 2005ab) to be included in the 
present compilation. Table lAll contains this new compila- 



tion and Appendixl^contains detail of its construction and 
references. It also provides a procedure to estimate stellar 
masses from the absolute K or V-band total magnitudes of 
a galaxy, given its Hubble Type. 

The local environmental galaxy density parameter (p; 
units of Mpc"^) from TuUy (1988a) is included in Table IXT] 
A smaller smoothing scale (cr = 0.5) than employed in 
the original catalogue is used to increase the density res- 
olution. Density p values ranging from 0-1, 1-5, 5-10, and 
> 10 Mpc""^ generally correspond to field, small group, large 
group (or small cluster), and cluster environments, respec- 
tively. 

Errors on the p values were approximated with the fol- 
lowing formula: Ap = \{Fp) — p\/p, where _F(m-M) is from 
TuUy (1988b) and is the expected number of galaxies over 
the number actually observable due to an increasing incom- 
pleteness fraction with distance. The average uncertainty of 
the present sample is Ap ~ 0.7. The TuUy (1988a) cata- 
logue relies on recession velocities for distances, therefore 
individual peculiar motions could influence p, especially in 
low-density environments where the low counts are prone to 
small number statistics. 

NGC 821 one of the most distant galaxies in the present 
compilation, thus the corresponding uncertainty in p (p = 
0.32 ± 2.0) is too large to prove its isolated nature. How- 
ever, confirmation of the isolated environ of NGC 821 was 
independently determined by Reda et al. (2004; and others 
prior, referenced in Reda et al.) from the larger LEDA galaxy 
catalogue (Paturel et al. 2003). Since the TuUy (1988a) 
NGC 821 p is likely under-estimated due to catalogue in- 
completeness, an upper-limit of p = 0.95 ± 0.5 was as- 
sumed. This value is the TuUy (1988a) density value of the 
nearby galaxy NGC 4594, which, according to an appropri- 
ate NED query, meets the isolated criteria of Reda et al. 
(2004). Other galaxies with distance moduli greater than 
m-M= 31.5 (NGC 1407, NGC 3268 and NGC 3258) were 
omitted from the density analysis as the p values are very 
likely to be underestimated. 

Many GC system properties are known to correlate with 
the mass of the host galaxy. It is therefore important to un- 
derstand selection biases relating to galaxy mass and envi- 
ronment. Indeed, Figure [T3l shows that the current sample is 
strongly biased towards very massive galaxies in the highest 
density environments. Part of this results from the preferen- 
tial examination of the most populous GC systems in order 
to compensate for the relatively large distances to galaxy 
clusters, but also reflects the fact that high-density regions 
are more likely to contain the most luminous galaxies. 

Before presenting the results, the use of the GC T- 
parameters (Zepf & Ashman 1993) instead of the galaxy 
luminosity- normalised GC specific frequency (Sjv) should be 
motivated. The T-parameter quantity, Tjv, is defined as the 
number of GCs normalised by its host galaxy stellar mass 
(see Tab.[7)|. The Tjv carries more physical meaning than Sjv 
because it is designed to account for stellar population dif- 
ferences (e.g. an old elliptical versus a young spiral galaxy) 
that bias optical photometry traditionally used for Sjv cal- 
culations (Zepf & Ashman 1993). 

To improve clarity in the following discussion, two new 
quantities are introduced as galaxy baryon and halo mass 
analogues to the Sjv and Tjv quantities. Table [7] contains the 
full definitions of the new (Ujv and Vjv) and old GC mass- 



© 0000 RAS, MNRAS 000.[llf23l 



12 L. Spitler et al. 





1 

XS) N821 


• 


11.5 


- • • ■ ■ . 1 • 1 

- ° ■ • ' . H 


m ' 

. \ m 

m 
















" 11 






10.5 







p(Mpc"^) 



Figure 13. Galaxy stellar mass versus local galaxy density from 
TuUy (1988a). Large symbols representing galaxies with GC sys- 
tems studied are as indicated in the figure legend. Small blue 
points and open red circles are nearby (m-M< 31.5) TuUy (1988a) 
catalogue spiral and elliptical galaxies, respectively. A compari- 
son between the two datasets suggests observational biases are 
evident in the present sample: GC systems of massive galaxies 
are preferentially studied in high-density environments. Galaxy 
stellar masses are from K-band photometry (see Appendix A for 
details). Errors on p reflect incompleteness fractions in the TuUy 
(1988a) catalogue (see text). 



Table 7. Globular cluster specific frequencies. Sjv is from Harris 
& van den Bergh (1981) and Tjv is from Zepf & Ashman (1993). 
Subscripts of N, blue, and red on the variables correspond to val- 
ues calculated with the entire GC system, blue GC subpopulation 
and red subpopulation numbers, respectively. 



Variable Definition Notes 

Sjv = NgcW°-^<-^''v + ^5) Mv is galaxy total V-band mag. 

Tjv = j^^'^igii ^* gE'l^'Xy stellar mass. 

Uiv = MtJ^Z'!^/W^ ^baryon is galaxy baryon mass. 

'^^ ^ Mhlu.fw^'^ A^haio is galaxy halo mass. 



normalised quantities. Table IA2I presents the GC number 
specific frequencies for the galaxies studied here. 

The present contribution to this data compilation, the 
isolated elliptical NGC 821, is shown as a ringed, magenta 
circle in subsequent figures. 

4.3 GC T-parameter Trends with Galaxy Stellar 
Mass and Environment 

Galaxy total V-band magnitudes are traditionally employed 
to approximate the stellar masses (Af, ) for deriving T- 
parameter values (e.g. Zepf & Ashman 1993; RZS05). How- 
ever, optical photometry is not a perfect proxy for M, . For 



instance, one correlation that should be accounted for when 
approximating M, from photometry is the galaxy mass- 
metallicity trend, which may account for up to 1/3 of the V- 
band T-parameter trends (e.g. RZS05). Infrared light from 
filters such as the K-band, are much less sensitive to metal- 
licity and therefore preferred for mitigating the galaxy mass- 
metallicity trend. Age corrections as well as effects from dust 
extinction are also much smaller in the K-band relative to 
V (e.g. Olsen et al. 2004). Another important advantage is 
that the K-band data come from a homogenous sample of 
the 2MASS (Jarrett et al. 2003) catalogue, rather than the 
popular, but heterogeneous compilation of V-band photom- 
etry from the RC3 (de Vaucouleurs et al. 1991). 

The text in Appendix A explains the adopted conver- 
sion between galaxy V and K-band photometry and stellar 
mass. Briefly, total stellar masses are the sum of the stellar 
mass in the bulge and disk components of a galaxy (ellip- 
ticals were assumed to have no disk). Appropriate mass-to- 
light ratios for assumed ages of the two components were 
taken from simple stellar population models (Bruzual & 
Chariot 2003) and bulge-to-total ratios for a given galaxy 
Hubble Type are from Graham & Worley (2007). 

GC subpopulation T-parameter values were computed 
using both V and K-band M, estimates and are presented 
for comparison in Figures [14] and 1151 respectively. Overall, 
the trends are very similar in the V and K-band versions, 
suggesting the assumptions made to estimate M« are valid 
for the current sample. Only a slight increase in apparent 
scatter in the K-band version of the Tuue versus log A/* plot 
for log M* < 11.2 is observed, which is likely just a small 
number effect. Also, the stellar mass estimates of the Scd 
spiral galaxy NGC 3556 demonstrate the advantage of using 
K-band photometry to estimate M,. This galaxy's V-band 
M* estimate (log A/* = 11.0) is 0.3 dex greater than the K- 
band estimate (log = 10.7), which is consistent with the 
idea that the V-band stellar mass is over-estimated due to 
a significant V-band emission from stars with ages < 5 Gyr. 
Thus, given such advantages, the K-band photometry will 
be used in the remaining figures and discussion regarding 
M.. 

4-3.1 Galaxy Stellar Mass Trends 

Aside from one obvious outlier (see the Fig. 1151 caption), 
it is clear that the larger sample considered here confirms 
RZS05, where a positive trend between Tbiun and host Af* 
was found. This trend implies that more massive galaxies 
exhibit larger stellar-mass normalised blue GC numbers. 

As pointed out in RZS05, this result is inconsistent with 
the major-merger model of Ashman & Zepf (1992), where 
the blue subpopulation of massive ellipticals is actually the 
sum of two blue GC subpopulations intrinsic to the two 
merging spiral galaxies. In this model, to reach the high 
Tbjue values of massive ellipticals, the merging spiral Tuue 
values must have similarly high values, contrary to what is 
confirmed here. Note, Burkert, Naab & Johansson (2007) 
have recently reported that their theoretical predictions for 
major galaxy-galaxy merger events are inconsistent with the 
observed central kinematic structure of massive elliptical 
galaxies. 

RZS05 speculate the positive Tuuc trend may relate to 
the timescales of GC formation in a simply biased galaxy for- 
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Figure 15. Same as Figure [141 except the T-parameters and M* were computed using galaxy K-band photometry. As detailed in the 
text, the p-statistics from Kendall tau correlation tests are given in the respective panels. Dotted lines are linear fits to the data. One 
outlier is evident and is labelled in the Tblue^^* plot: the isolated spiral NGC 7331 (Rhode et al. 2007). 



mation model (Santos 2003). Under such a model, massive 
galaxy potentials hosted the earliest sites of GC formation, 
thus more massive galaxies may have had more time to pro- 
duce blue GCs (see discussion in RZS05; Santos 2003; also 
Brodie & Strader 2006). 

Noticeably absent from the RZS05 Tuue study (due to a 
scarcity of good total GC population estimates at the time) 
were massive central cluster galaxies (e.g. Virgo Cluster M87 
and Fornax Cluster NGC 1399), which have been shown 
to host enhanced galaxy luminosity-normalised GC num- 



bers relative to other galaxies (e.g. Harris & van den Bergh 
1981). Included here are four galaxies located at the centre of 
the three nearest galaxy clusters: Virgo, Fornax and Antlia 
(NGC 3258 and NGC 3268 are each located in a central 
sub-group within the Antlia cluster). These galaxies show 
the largest T-parameter values of the sample (Tuue > 10 
and Tred > 6), confirming previous observations. The log- 
arithmic scale of the T-parameter values in the plot helps 
support the idea that these GC systems are not actually 
unique, but simply at the extreme end of a general trend. 
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Figure 14. V-band blue and red GC T-parameters versus host 
galaxy stellar mass left panels) and local galaxy number 

density (p; right panels). Tblue ^.nd Tred E're the number in the 
respective GC subpopulation normalised by its host galaxy stel- 
lar mass. Mt is computed using total V-band galaxy luminosity. 
A trend for increasing Tuue with galaxy M« is apparent, while 
the Tred plot shows a weaker trend. Galaxies with high p values 
generally show larger T-parameter values. 

The Tred data shows evidence for an positive analogous 
trend with M*, but with more scatter. RZS05 mentioned, 
but did not discuss the Tred trend. Again, this appears to 
imply that massive galaxies formed relatively more red GCs 
compared to lower mass galaxies. Larger scatter may indi- 
cate the formation process of red GCs was subject to more 
stochastic variations than the processes responsible for blue 
GC subpopulations. In contrast to the early-type galaxies, 
the Tred values for spiral galaxies show no significant corre- 
lation with galaxy mass or p. This may suggest the red GC 
subpopulations of spirals were formed in a different man- 
ner than those of early-type galaxies (e.g. effects of secular 
evolution vs. merger-driven spheroid formation). 

As mentioned previously, the subpopulation T- 
parameter values of NGC 821 are fairly consistent with those 
of similar-mass spiral and elliptical galaxies. 

Least-square linear fits were performed on the data 
(clipping appropriate outliers) and are presented in Fig- 
ure [TS] and Table [S] Non-parametric, Kendall tau rank cor- 
relation coefficients were calculated for the data and their 
associated p-statistics are presented in the respective panels 
of Figure [T^ The p-statistic is the probability of the null hy- 
pothesis that the data are not correlated. According to the 
Kendall measure, both M, trends show a significant trend 
at a J5 99.99% confidence level. 

4-3.2 Local Galaxy Environment Trends 

On the right panels of Figure 1151 the GC subpopulation 
T-parameters are plotted against the local galaxy number 



Table 8. Linear fits to the T-parameters versus stellar mass and 
local density data presented in Figure [TSl The p-statistic is from 
a non-parametric Kendall tau test for the significant of a linear 
correlation. 



GC 


log T=a 


log M. + (3 


Kendall tau 


subpop. 


a 


/9 


p-stat. 


Blue 


0.80 ±0.17 


-8.3 ± 1.9 


0.0006 


Red 


0.66 ±0.22 


-7.0 ±2.5 


0.0102 


GC 


log T=a log p + /3 


Kendall tau 


subpop. 


a 


/9 


p-stat. 


Blue 


0.32 ±0.17 


0.29 ± 0.14 


0.2387 


Red 


0.27 ±0.19 


0.15 ± 0.15 


0.1939 



density, p, which is a proxy for the local, underlying mass- 
density. While p is not the only or even the most sensitive 
parameter for gauging environmental infiuence, it is useful 
for nearby galaxies and indeed provides a quantitative mea- 
sure of local environmental density (e.g. Blanton & Berlind 
2007). 

Compared to the T-parameter plots against mass, the p 
plots in Figure [T5] only hint at a positive trend. Supporting 
this qualitative assessment of the data, a Kendal tau test 
indicates p correlations at the ~ 99.8% confidence level and 
the linear fits in Table [S] are only significant at the 1 — 2a 
level. Nevertheless, on average, the T-parameters are lower 
in low-density environments and larger for higher values of 
p, with galaxies in a given environment spanning a factor of 
10 in T-parameter values. The isolated elliptical NGC 821 
is among the lowest T-parameter values in the sample. 

It is perhaps reasonable to ask if either M« or p rep- 
resents the dominant physical quantity eflecting the T- 
parameter values. As massive galaxies naturally tend to be 
found in higher-density environments, it is unlikely that 
completely separating galaxy mass and environmental ef- 
fects on GC system properties is possible or even an appro- 
priate course of action. Nevertheless, we can consider galax- 
ies of a similar mass or environment and ask whether or not 
the observed scatter in GC numbers is correlating with the 
third parameter. The present sample bias towards high M» 
galaxies in dense environments (see Fig. I13|) unfortunately 
makes any such inference very uncertain without a more 
complete dataset. 

Despite this caveat, each panel in Figure [16] shows the 
residuals of linear fits to the data in the corresponding panels 
of Figure 1151 For instance, the residuals fit to the Ttiue- 
Af» data in the upper left panel of Figure 1151 are presented 
in the corresponding panel of Figure [16] against the third 
parameter, p. A very shallow trend is visible in the right 
panels of Figure 1161 perhaps suggesting that a M* trend 
remains after environmental effects are removed. This seems 
to suggest a stronger dependence on M*. However, this may 
only refiect the poor, 2a fit of the T-parameter and p data. 

Thus tentative evidence is found that the T-parameters 
depend more on M« than p, although separating the two ef- 
fects is indeed challenging at the present. It is noted that the 
isolated elliptical NGC 821 shows a consistent GC number to 
other similar-mass galaxies in higher density environments 
(see Figure [TS]). Other measures of local environmental den- 
sity (e.g. group luminosity, group-centric distances) might 
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Figure 16. Residuals for the Tmub {top panels) and T^^^ {bot- 
tom panels) trends versus galaxy stellar mass (M*; right panels) 
and local galaxy density (p; left panels). Residuals are from lin- 
ear fits presented in Figure 1151 Residual RMS values are given. 
Each panel corresponds to the same panel of Figure 1151 with the 
third parameter not used in the fit presented in this figure on 
the X-axes. Only a hint of a trend is visible in the right pan- 
els, which would indicate M» represents the dominant physical 
property effecting T-parameter values. 
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Figure 17. Similar to Figure [141 except the U-parameter is plot- 
ted against stellar plus X-ray hot gas mass {Mi,aryon)- The U- 
parameter trends show only a slightly shallower gradient com- 
pared to the T-parameter data, thus the overall effect of including 
hot gas is minor. 



Table 9. Like Table [S] except the fits are to the Mi^aryon U- 
parameters, presented in Fig. 1151 



help distinguish whether host galaxy mass or local environ- 
ment has a stronger effect on GC system properties. Further 
discussion of the environmental trends is provided in H4.8I 



4.4 Galaxy Baryonic Mass Trends 

Before interpreting the T-parameter trends in terms of GC 
formation efficiency, the normalising quantity in the T- 
parameters, M« , must be understood physically. Total lumi- 
nosities of galaxies bare a direct connection with the stellar 
mass content of a galaxy, once the stellar populations are 
understood. However, processes such as strangulation can 
lead to a decreased field star formation efficiency in cluster 
galax;ies while having no effect on GC numbers. This effect 
can induce or strengthen a positive T-parameter trend with 
host stellar mass, thus normalising GC numbers with a total 
haryon mass and computing Ujv (see Tab. [7)l may be more 
relevant if a connection between GC and galaxy formation 
is to be attempted. 

McLaughlin (1999) explored this idea and found evi- 
dence for a universal fraction of a galaxy's baryon in its 
GC system (which would roughly correspond to constant 
U(,iue-|-Ur-ed values) , across a range of galaxy luminosities. 
Since the McLaughlin (1999) study, many of the GC system 
number estimates have decreased (e.g. Rhode et al. 2007), 
thus an updated analysis is pertinent and is presented in 
Figure [17] and Table [9] U-parameters in this figure use a 
galaxy mass estimate that includes both stellar masses from 
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K-band photometry and the X-ray hot gas mass fractions 
from McLaughlin (1999), which is reproduced here: 

M,„, = 0.55M,(^^^)i'^ (2) 

where Ly is the V-band luminosity of the host galaxy. As 
galaxies show a fairly large intrinsic scatter in their X-ray lu- 
minosities at a given stellar mass (e.g. Mathews et al. 2006), 
the subsequent baryon mass analysis based upon Equation 2 
is only useful for studying overall trends. 

Interestingly, the U-parameter trends in Figure [T7] per- 
sist with only a slight decrease in the T-parameters gradients 
in Figure [151 The slightly shallower trends are caused by a 
greater hot gas fraction in massive galaxies. Linear fits in 
Figure [TTl and Table [9] support this qualitative assessment. 

The universal GC mass {Mac) fraction of the host 
baryon mass {M},aryon) found by McLaughlin (1999) was 
= 0.0026, where e'' = Mac /Miaryon- The bottom panel 
of Figure [TSl shows the total mass of GCs against the host 
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Figure 18. GC mass (Mqc) against host galaxy baryon mass 
i^baryon)- Red, blue and all GCs of a given galaxy are presented 
in the top, middle and bottom panels, respectively. The dotted 
lines correspond to constant values of = Mqc /^baryon = 0-01 
and = 0.0001. Solid lines are least-squares fits to the GC data and 
the dotted line in the lower panel corresponds to the "universal" 
e'' = 0.0026 determined by McLaughlin (1999). Newer GC num- 
ber estimates have decreased due to better account of contamina- 
tion, thus the data presented here is inconsistent with a univer- 
sal GC to galaxy baryon mass fraction as found by McLaughlin 
(1999). 



baryon mass, where the GC system masses are the product 
of the GC population number and the average CC mass of 
2.4 X IO^Mq used in McLaughhn (1999). Included in the 
Figure are lines of constant e*" and the McLaughlin (1999) 
universal value of e** = 0.0026. 

In the present sample, e*" of massive galaxies are con- 
sistent with the universal McLaughlin (1999) value. How- 
ever, for lower mass galaxies, the GC masses clearly show 
lower values than the "universal" e*". A least-squares fit to 
the data yields the equation log Mac = 1-59 ± 0.16 log 
Mbaryon — 9.5 ± 1.8, which is inconsistent with a constant 
(i.e. slope of 1) at the ~ 4(t level. Both GC subpopulations 
show similar linear trends with only a zeropoint offset. 

To reach a constant value of e'', a factor of 5 more 
baryons would need to be added to the mass budget of mas- 
sive ellipticals. Otherwise, the baryonic masses of low-mass 
galaxies could be over-estimated by a similar factor. Both 
extremes seem unlikely. As this analysis does not include 
dust or gas too cold for detection in X-rays, it is techni- 
cally incomplete and a more perfect analysis would involve 
a detailed case-by-case accounting of an individual galaxy's 
baryon content. However, these contributions should be rel- 
atively minor for massive galaxies and an increase in the 
baryon content of low-mass galaxies would only strengthen 
the U-parameter trends observed here. 

This result can be corroborated with the use of dynam- 
ical masses (Mdyn)- Galaxy Mdyn estimates are from the 
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Figure 19. Comparison between dynamical and total baryon 
galaxy masses, whose derivation is discussed in the text. Com- 
parison to the 1:1 line indicates the two mass estimators yield 
consistent results. 

Cappellari et al. (2006) form of the viral theorem is used, 
with central velocity dispersions are from LEDA (Paturel et 
al. 2003) and K-band effective radii from 2MASS 20 mag 
arcsec"^ isophotal radii, converted to effective radii with 
an empirical relationship from data in Jarrett et al. (2003). 
Figure [19] shows that the dynamical masses are compara- 
ble to Mtaryon, Suggesting that the baryon masses derived 
from photometry are a good indicator of the total baryonic 
content. 

Hence, contrary to the conclusions of McLaughlin 
(1999), this analysis suggests GC formation efficiency in- 
creases with the host baryon mass. The main cause for the 
new interpretation is that new GC number estimates have 
become available, which were derived using improved con- 
tamination estimates and better constrained radial profiles 
of the GC systems (see e.g. Rhode et al. 2007). 

4.5 Galaxy Total Halo Mass Trends 

A final quantity worth exploring for the GC number normali- 
sation mass is galaxy total halo mass, Mhaio- Mandelbaum et 
al. (2006) performed a galaxy-galaxy weak lensing analysis 
of SDSS galaxies and derived mean halo masses for galaxies 
spanning a range of stellar masses (log M« = 9.9 to = 11.6). 
They found a non-linear, positive gradient between galaxy 
stellar and halo mass that steepens significantly for galax- 
ies with log M* > 11. Although such pioneering work must 
be used with caution, the Mandelbaum relationship can be 
used to derive total halo masses and hence new normalised 
GC numbers, the V-parameters (see Tab.[7|, calculated. 

A few caveats should be considered first. As the lens- 
ing sample of Mandelbaum et al. (2006) occupy a range of 
redshifts (0.02 < z < 0.35), the halo masses might be un- 
derestimated for the local galaxy sample, because halos (and 
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Figure 20. GC halo mass normalised numbers (Vjv) against galaxy halo mass (Mhalo), which is derived from the weak lensing analysis 
of Mandelbaum et al. (2006). V-parameter values are normalised to a halo mass of log M^aio = H a^nd show no significant correlation 
with either Mf^^lo or p. This suggests GCs formed in direct proportion to the total halo mass of a galaxy. See text for detailed discussion. 



their GC systems) may grow with time. However, Parker et 
al. (2007) find little or no evidence for halo mass growth over 
the redshift range 0.2 < z < 0.8, and if this trend can be 
extended to the present epoch, would suggest the Mandel- 
baum et al. (2006) halo masses are roughly valid for nearby 
galaxies. Another potential concern is that halos likely show 
intrinsic variations for a given galaxy stellar mass, thus the 
halo masses assumed may not be accurate representations 
of the specific galaxies considered here. Mandelbaum et al. 
(2006) must stack their sample of galaxies to detect a lens- 
ing signal. Their work therefore gives a good average halo 



in a statistical sense, but will not apply to an individual 
galaxy. Mathews et al. (2006) provide some evidence this 
variation is fairly significant for a given galaxy M* , thus the 
subsequent analysis is only valid for overall trends. As to the 
robustness of the Mandelbaum et al. halo masses, it is noted 
that other studies have produced similar results. Parker et 
al. (2007) confirm the Mandelbaum et al. (2006) halo masses 
for L* galaxies with a deeper weak-lensing analysis. Halo oc- 
cupation statistical modelling by van den Bosch et al. (2007) 
predicts a similar scaling of halo and central stellar mass. 
The V-parameter values are arbitrarily normalised to 
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a halo mass of log Mhaio ~ 11 (see Tab. [7]) and are pre- 
sented in Figure [20] Mandelbaum et al. (2006) relationships 
for early and late-type galaxies were applied to the appropri- 
ate galaxy type for galaxies log M« > 11. Below this stellar 
mass, Mandelbaum et al. find no significant difference be- 
tween the halo masses of early and late-type galaxies, and 
so the early-type relation was used in this range. Figure [20l 
shows a striking contrast to the T and U-parameter plots 
discussed previously: the positive trend with galaxy mass 
is gone, suggesting GCs (of both subpopulations) formed 
in a constant proportion to the host halo mass, irrespec- 
tive of environment or host galaxy type. The Mhaio data 
are presented in alternate form in Figure [21] and are sta- 
tistically consistent (according to least-squares fits) with 
a constant GC mass fraction of the host halo mass. The 
best-fitted line to the total GC system properties is: log 
Mac = 0.94 ± 0.06 - 3.77 ± 0.76. 

The idea that GCs formed in direct proportion to the 
mass of its host halo is attractive for a number of reasons (see 
discussion in Blakeslee, Tonry & Metzger 1997 and Blakeslee 
1999). Such a model gels with observational evidence that 
most GCs formed in the very early (> 10 Gyr) Universe, 
conceivably before galaxy star formation histories began to 
be influenced by processes like AGN feedback or strangu- 
lation. In this framework, central massive galaxies such as 
NGC 4486 or NGC 1399, show "enhanced" GC populations 
(see i ]4.3.1|l because they are associated with an extremely 
massive dark matter halos. A useful consequence is that a 
direct relationship between GC system numbers and total 
halo mass could be derived. Such a relationship would help 
characterise the undetected or removed baryon content of a 
galaxy. 

A few theoretical predictions exist for this analysis 
(Santos 2003; Diemand et al. 2005; Kravtsov & Gnedin 
2005). In i |4. 71 below, the blue GC results will be compared 
to the Diemand et al. (2005) predictions for a variable reion- 
isation epoch. Santos (2003; see also RZS05 and Rhode et al. 
2007) predicted that galaxy "biasing" will naturally lead to 
greater blue GC numbers in more massive galaxies (contrary 
to what is found here) assuming blue GC formation is uni- 
versal truncated by a mechanism like reionisation. Santos 
(2003) suggested the sites of massive galaxies were biased 
to host the first GC formation regions, thus these galaxies 
would have had a longer opportunity to produce blue GCs 
before a universal truncation. 

Kravtsov & Gnedin (2005) reported a correlation be- 
tween the mass of a GC system and its host halo mass for 
simulated halos between log Mhaio = 10.8 and = 11.4. Al- 
though their analysis is not valid for the halo mass range 
considered here. Figure [21] shows the correlation extrapo- 
lated from lower masses. Kravtsov & Gnedin (2005) point 
out that the Milky Way (here with log Mhaio ~ 12.5) falls on 
their extrapolated relation. Supporting this observation, the 
present analysis suggests all galaxies in the halo mass range 
log Mhaio 12 to 13 are roughly consistent with their extrap- 
olated predictions. Only the most massive galaxies appear 
to deviate. 




z" ; I I I i -r'T"" i I I I I I I I I I I I I 




12 12.5 13 13.5 14 

log(M,„J (MJ 

Figure 21. Total mass in GCs (Mqq) against host galaxy halo 
mass (Mhaio)- The layout is as in Figure [TSl with data from Fig- 
ure 1201 Dotted lines correspond to arbitrary constant values of 
Mac /Mhaio = e'^ = 1 X 10""' and = 1 X 10~^. The data are 
consistent with constant GC mass fractions. Mean mass fractions 
of blue, red and all GCs are e'' = 2.0 X 10-^,1.3 X 10"^, and 
3.2 X 10~^, respectively. Dashed line is the extrapolated predic- 
tion from theoretical analysis of Kravtsov &; Gnedin (2005) . 

4.6 Comparison between Galaxy Halo and 
Baryon Mass Trends 

Although galaxy stellar to halo mass fractions [M, / Mhaio) 
vary with host galaxy mass (e.g. Mandelbaum et al. 2006), 
on scales of galaxies the baryon fraction might roughly re- 
semble the universal baryon fraction (i.e. Mtaryon/Mhaio ~ 
Qt/^n, ~ 0.17; WMAP3, Spergel et al. 2007). If a universal 
baryon fraction is applied to the Mcc-Mbaryon data in Fig- 
ure [18] to make a Moc-Mhaio plot analogous to Figure [21] 
only an X-axis shift to larger halo masses would occur and 
the form of the Mac / Mtaryon trend will directly translate 
into to a Mac /Mhaio trend. This apparently contradicts the 
roughly constant Mac /Mhaio shown already in Figure [51] 
A possible explanation for this "discrepancy" is that the lo- 
cal baryon fraction decreases with increasing host halo mass. 
Although this idea is interesting, a few analytical checks and 
further discussion are needed before making any conclusions. 

It is first worth reviewing the empirical relationships 
used to derive the Mbaryon and Mhaio values. Total galaxy 
baryon masses are from the sum of total stellar mass (from 
K-band photometry; see plus an X-ray hot gas compo- 
nent (between 40% and 0% of M, for high and low-mass 
galaxies, respectively; see Eq. 2 in MA\ and McLaughlin 
1999). The hot gas relation is normalised to NGC 4486 
(M87), for which McLaughlin (1999) estimate ~ 40% of the 
baryons are in hot gas within 100 kpc of the galaxy centre. 
To make Mac /Mbaryon ~ constant, this percentage would 
have to increase to an unlikely value of ~ 90% (i.e. thereby 
increasing baryons by a factor of five, see ^4.4|l . suggesting 
the current baryon calculations are reasonable. It was also 
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found that the Mtaryon are consistent with mass estimates 
inferred from central galaxy dynamics (see Fig. I19|l . 

The contradiction may also be alleviated if the Man- 
delbaum et al. (2006) halo masses are under-estimated for 
low-mass galaxies or over-estimated for high-mass galaxies. 
Interestingly, Mandelbaum et al. (2006) find that for log 
M* < 11, early and late- type galaxies (classified as such by 
their light profiles) have similar Alhaio values, on average, 
but for log A/, > 11 early-type galaxies are hosted by sig- 
nificantly more massive halos than those surrounding late- 
type galaxies. Mandelbaum et al. point out that massive 
early-type galaxies tend to be found at the centres of galaxy 
groups or clusters and speculate the lensing masses of mas- 
sive early- type galaxies actually refiect group/cluster mass 
potentials rather than galaxy-sized halos. If this interpreta- 
tion of Mandelbaum et al. is correct, then the M^aio values 
of massive early- type galaxies in Figure [20] and [21] may not 
be comparable to those galaxies of log A/, < 11. Applying 
the late-type trend (i.e. from non-central galaxies) trend to 
all galaxies of log A^« > 11 increases their halo V-parameter 
values by an amount that may remove the discrepancy with 
the Mbaryon U-parameter trends. This seems to support the 
Mandelbaum et al. interpretation of the observed halo mass 
difference between early and late-type galaxies. 

On the other hand, four of the galaxies considered here 
are actually located at the centres of massive groups or 
galaxy clusters. It is also very likely that very massive, non- 
central eUipticals (e.g. NGC 4472 and NGC 3268) must have 
once been at the centre of a large galaxy group, which has 
since been incorporated into a galaxy cluster. As these mas- 
sive ellipticals host enhanced GC systems for their bary- 
onic masses (see i^4.4|) . it is worth asking whether there ex- 
ists a connection between the GC system and the larger 
cluster-sized halos they are (or once were) at the centres of. 
Blakeslee et al. (1997; see also Blakeslee 1999) favoured such 
a connection after they found that the Sjv values of massive 
elliptical galaxies at the centres of clusters correlate strongly 
with the host cluster mass. 

To remove the apparent discrepancy between the Uiv 
and Vat values, the Uiv values of the massive galaxies must 
be brought down to the observed levels for "normal" low- 
mass galaxies, thus more baryons need to be added to the 
central galaxy's mass budgets. Assuming massive elliptical 
GC systems are or once were associated with a cluster-sized 
halo, a more appropriate M.t,aryon estimate might include 
baryons associated with the intra- cluster medium. The con- 
tribution from an intra-cluster baryon component for a cen- 
tral cluster galaxy mass budget can be crudely estimated 
here. Gonzalez, Zaritsky & Zabludoff (2007) found that the 
luminosity of the intra-cluster light is roughly greater than 
5 times the light in the brightest cluster galaxy, so includ- 
ing this component in the cluster baryon budget would add 
roughly a factor of ~ 5 or more baryons. This factor is very 
similar to the increase in baryon mass needed in the most 
massive ellipticals to produce a constant value of Uat for all 
galaxies (see i]4.4p . Baryons from intra-cluster gas will push 
this factor beyond five, which might be counter-balanced by 
allowing for a GC mass component from intra-cluster GCs. 
Though a future, more detailed investigation is needed, this 
analysis lends some support to the idea that central galaxy 
GC system numbers are indeed proportional to total cluster 
mass as found in Blakeslee et al. (1997) and Blakeslee (1999). 



Furthermore, the inclusion of baryons from the intra-cluster 
medium surrounding central galaxies may remove some, if 
not all, the discrepancy between the Uat and Yn trends (im- 
plying the baryon fraction, on the scales of galaxy clusters, 
is equivalent to the universal baryon fraction). This again 
suggests that GCs formed in proportion to the total halo 
mass available, be it from a cluster or galaxy-sized halo. It 
is therefore appropriate to use the Mandelbaum et al. (2006) 
early-type (i.e. central) relationship on these massive ellip- 
tical galaxies, as in Figures [20l and [2T] 

Clearly further work is needed on many fronts. GC sys- 
tem numbers show at least some promise as indicators of 
halo mass, although individual halo estimates of galaxies in 
the local Universe (e.g. from GC dynamics, Romanowsky 
2006) are needed to confirm and better quantify such a cor- 
relation. Carefully determining the baryon content of indi- 
vidual galaxies and comparing their GC system numbers to 
other galaxies will help determine the extent of the possible 
GC number and baryon/halo mass discrepancy and whether 
the local baryon fraction varies with host galaxy mass. GC 
systems of dwarf galaxies are very interesting systems to add 
to this analysis because they are thought to be dominated 
by dark matter, thus they might host many more GCs than 
otherwise would be predicted from their stellar masses. 



4.7 Constraining Reionisation 

Although contrary theories exist (e.g. Gen 2001), several 
workers have speculated how GC metallicity bimodality can 
result from the truncation of GC formation during cosmic 
reionisation (e.g. Santos 2003; Beasley et al. 2002; Moore et 
al. 2006). In this picture, metal-poor GC formation halted 
when reionisation increased the ambient temperature be- 
yond the threshold for efficient star formation. Some period 
of time later, after gas had become enriched with metals and 
cooled, a second GC formation epoch commenced resulting 
in the second GC metallicity subpopulation. Though the de- 
tailed physics of GC formation is not well-constrained, GCs 
provide ideal tracers of dark matter from very early stages 
of galaxy formation (Diemand et al. 2005), thus they could 
provide observational constraints on the nature of reionisa- 
tion (Santos 2003; Diemand et al. 2005; Moore et al. 2006; 
Bekki 2005 and for intra-cluster GCs see Bekki & Yahagi 
2006 and Power, Drinkwater & Silk 2007). 

Diemand et al. (2005) predicted that galaxies with the 
largest Vuue values will tend to reside in regions where reion- 
isation occurred later: i.e. locations where metal-poor GC 
formation was more prolonged. In the context of this model, 
the constant Vuue values for all galaxies (see H4.5P suggests 
the epoch was roughly coeval throughout the Universe. A 
fiat trend with environmental density may also imply that 
reionisation was not significantly biased to start in a par- 
ticular type of environment (e.g. Weinmann et al. 2007). 
Detailed modelling of the range of galaxy halo masses ex- 
amined here is required before the observations can provide 
tighter constraints on the exact epoch of reionisation and 
the nature of its propagation throughout the Universe. 
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4.8 Globular Cluster Systems and their 
Environment 

Brodie & Strader (2006) discussed recent observational evi- 
dence that may imply GCs formed in the proto-galactic frag- 
ments before galaxy assembly began. Under such a scenario, 
GC numbers may provide a measurable proxy for the initial 
conditions of a given galaxy (e.g. the total available mass or 
environmental conditions) . Thus GC numbers could serve as 
a basis quantity for comparison between different galaxies to 
isolate evolutionary processes that shape galaxies into their 
present-day form. Before GC systems can be used as such, 
the fundamental relationship between their own properties 
and the initial conditions (e.g. available mass and their envi- 
ronment at formation), must be understood independent of 
their host galaxies. The objective of this section is to deter- 
mine whether or not GC system numbers (i.e. GC formation 
efficiencies) show intrinsic differences with the local environ- 
mental density, independent of GC system number trends 
with the host galaxy mass. Galaxy baryon masses are pre- 
ferred over M, for the normalising quantity in the GC spe- 
cific frequencies, because Mtaryon (like GC numbers) is not 
influenced by certain environmentally-dependent processes 
that will effect a galaxy's total stellar mass content. 

In the past, the answer to the question whether or 
not GC system formation is sensitive to the local environ- 
ment, has tended towards the affirmative: GC systems were 
thought to be intrinsically more populous in higher density 
environments, implying GC formation efficiency was some- 
how higher in these environments (e.g. Ifarris 1991; West 
1993; Kumai, Hashi & Fujimoto 1993). Analysis presented 
in Section 14.3.21 largely confirms Kumai et al. (1993) and 
West (1993), who found a positive trend between Sjv and 
the TuUy (1998a) environmental density parameter. How- 
ever, the robustness of this result is questioned here. Due to 
a bias in the early-type GC system number estimates avail- 
able, both the present and previous studies (see Fig. 4 in 
Kumai et al. 1993; a similar bias exists in the West 1993 
sample) do not have access to GC systems of low-mass clus- 
ter galaxies, thus there is not sufficient leverage to separate 
mass from environmental effects on the GC system numbers 
(see gXlJ. 

Such studies have therefore not presented solid evidence 
for a Sat (or Tjv and Uat; see tj4.4|l trend with local environ- 
mental density once trends with galaxy mass are removed. 
Nor has contrary evidence been presented, thus the extent of 
an environmental impact on non-central GC systems (inde- 
pendent to the host mass) is subject to debate. However, the 
tentative result for a constant halo mass-normalised number 
of GCs (i.e. Vjv) over the entire range of galaxies masses con- 
sidered here, is consistent with a model where GC system 
numbers show no or very little environmental dependence. 

Among giant ellipticals located at the centres of galaxy 
clusters, Blakeslee et al. (1997; see also West et al. 1995; 
Harris, Harris & McLaughlin 1998) found some evidence for 
an increase in Sjv with increasing number of nearby galaxies, 
i.e. environmental density. In addition, they found a strong 
correlation between central galaxy's Sjv values and the mass 
of the host cluster. It is noted that there exists an empirical 
trend where more massive clusters tend to show larger cen- 
tral galaxy number density values (e.g. Bahcall 1981). As the 
central galaxies studied by Blakeslee et al. (1997) occupy a 



very narrow range of stellar masses, a Sjv trend with mass or 
environment can be caused by a variation in the GC and/or 
field star formation efficiencies. Blakeslee et al. (1997; see 
also Blakeslee 1997) and others (Blakeslee 1999; Harris et 
al. 1998; McLaughlin 1999) adopt the view that the cen- 
tral galaxies of more massive clusters become increasingly 
under- luminous (from an increased inefficiency of field star 
production) and therefore central GC system numbers scale 
with the total cluster halo mass and/or local environmental 
density. It is again apparent that if environmental effects are 
to be isolated, baryon masses should be used instead of stel- 
lar masses, especially for massive ellipticals where hot gas 
contributes significantly to yibaryon- 

McLaughlin (1999) performed a careful baryon mass 
accounting for the central Virgo (M87) and Fornax 
(NGC 1399) cluster galaxies and found a very similar frac- 
tion of the galaxy's baryon mass is contained in these GC 
systems. Given that Virgo and Fornax clusters exhibit no- 
ticeably different masses and central galaxy number den- 
sities (see e.g. Jordan et al. 2007a), the McLaughlin (1999) 
results provide strong evidence that GC formation efficiency 
in central ellipticals is remarkably similar. Thus, while Sjv 
of massive ellipticals clearly correlates with the cluster mass 
and possibly with local environmental density (e.g. Blakeslee 
et al. 1997), it isn't clear that Ujv does. Furthermore, ob- 
servational evidence suggests that for massive ellipticals Vjv 
is roughly constant (Blakeslee et al. 1997; Blakeslee 1999; 
^43}. If the above were confirmed, GC formation efficiency 
is likely to be insensitive to local environmental density and 
hence GC system numbers can serve as a good indicator of 
the initial mass available to a host galaxy. 



5 SUMMARY AND CONCLUSIONS 

The present work on NGC 821 gives the best constraints 
on an elliptical galaxy's GC system properties located in a 
truly isolated environment. All measured properties of its 
GC system indicate a fair resemblance to similar-mass el- 
liptical galaxies in higher density environments. Its isolated 
condition has therefore produced no distinguishing features 
compared to ellipticals in cluster environments. 

This result is supported by a comprehensive compari- 
son of the GC systems in 25 galaxies spanning a range of 
Hubble Types, masses and environments. Two new GC spe- 
cific frequencies are introduced: Uat and Vat, which refer 
to GC numbers normalised by the host galaxy baryon and 
halo masses, respectively (see Tab. [T]). With better galaxy 
mass estimates using 2MASS K-band galaxy photometry 
and the inclusion of an hot X-ray gas component, the new 
compilation suggests that the galaxy baryon mass {Mbaryon) 
normalised number of metal-poor and metal-rich GCs (i.e. 
the GC subpopulation U-parameters: \5biue and Ured) in- 
creases with host galaxy baryon mass. This result is at odds 
with idea that GC formation is directly proportional to the 
Mtaryon of the host galaxy, as concluded in the study of 
McLaughlin (1999). The different conclusion reached here is 
likely due to improved contamination analysis used in recent 
GC system studies. 

Spiral galaxies in the sample generally show lower Tuue 
values compared to massive ellipticals, ruling out the pos- 
sibility that massive elliptical GC systems were formed via 
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the major mergers of spiral galaxies (Ashman & Zepf 1992) . 
Note, the T-parameter values of spirals are similar to low- 
mass ellipticals. Also, the spiral GC system studies used here 
are biased towards field spirals, and do not include spirals 
located on the outskirts of galaxy clusters. 

New weak lensing halo mass [Mhaio] estimates from 
Mandelbaum et al. (2006) were employed to investigate GC 
formation efficiency relative to the host halo mass. The Vjv 
data appears to be consistent with a constant GC number 
per unit halo mass, as found among massive central clus- 
ter galaxies by Blakeslee et al. (1997). GC system numbers 
might therefore provide an estimate of the total halo mass 
of a hosting galaxy. This analysis is compared to the U- 
parameter trends. There is a possible discrepancy between 
the two if the local baryon fraction does not decrease with 
increasing galaxy mass. Part of this problem may be re- 
moved if the GC system numbers of very massive ellipticals 
are directly proportional to the total halo mass of the cluster 
these giant ellipticals tend to be at the centres of. Caveats 
and implications are discussed in tj4.5l and >i4.6l 

The T and U-parameters are also found to correlate 
with local galaxy environment, as traced by a modified ver- 
sion of the TuUy (1998a) p parameter. Due to a bias in the 
present and past (Kumai et al. 1993; West 1993) galaxy sam- 
ples, it is currently difficult to determine whether Mtaryon 
or p has a greater inffuence on GC numbers. In central 
cluster galaxies, while Sjv may correlate with environment 
(e.g. Blakeslee et al. 1997), Uat likely does not (McLaughlin 
1999), supporting a model where GC formation efficiency is 
insensitive to the local environment. This model is further 
supported by the observation that Vjv values are roughly 
constant among massive, central ellipticals (Blakeslee et al. 
1997) and the assortment of galaxies studied here ( ^14.5^ . 

Theoretical work by Diemand et al. (2005) provides an 
interesting connection between the number of blue GCs and 
the local epoch of reionisation. They predict the relative 
number of blue GCs should be higher in regions of the 
Universe where reionisation took place later because the 
time available for blue GC formation was longer in duration. 
The GC systems considered in the present work may show 
a constant number of blue GCs per unit halo mass, Vuue, 
suggesting the reionisation did not vary significantly with 
host halo mass or the local environmental density. 
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APPENDIX A: DATA COMPILATION 

Literature compilations of robust GC system properties and 
their references are presented in Table \KT\ GC specific fre- 
quencies are tabulated in Table IA2I Galaxy Hubble Types 
were taken from NED and the local density parameter, 
p, was computed from the TuUy (1988a) catalogue, as de- 
scribed in i]4.2l Most distances were taken from the NED1E0 
database. Errors on T-parameter values are from published 
number uncertainties. If such numbers were unavailable, 
these errors reflect Poisson uncertainties of the GC num- 
bers. References for the distances and GC system properties 
are given in the Table \KT\ caption. 

To derive galaxy stellar masses from photometry, a 
mass-to-light ratio that varies with Hubble Type was used. 
Galaxy discs show younger stellar ages compared to bulge 
components, thus the total measured photometry was di- 
vided into these two components with empirical K and B- 
band (B-band is used here for the V-band galaxy photom- 
etry) bulge-to-total ratios from Graham & Worley (2007). 
Bruzual & Chariot (2003) theoretical mass-to-light ratios 
(assuming a metallicity of [Fe/H] — -1-0.09) for a young (5 
Gyr) and old (12 Gyr) stellar population were employed 
for the disc and bulge components, respectively. Elliptical 
galaxies have no disc component, thus their light is assumed 
to be dominated by an old stellar population. These mass 
estimates only include stellar mass and not the estimated 
gas produced during stellar evolution. Mass-to-light ratios 
for a given Hubble Type are presented in Table IA3I Note 
the values in Table IXSl are very similar to those used in the 
original definition of Tjv by Zepf & Ashman 1993, suggest- 
ing their empirical mass-to-light ratios vary due to stellar 
population effects. 

As discussed in H4.3I approximating galaxy masses with 
K-band photometry is preferred over V-band. Stellar masses 
in Table are from K-band photometry, which is from the 
2MASS Extended Source Catalogue (Jarrett et al. 2000) and 
the 2MASS Large Galaxy Atlas (Jarrett et al. 2003) when 
appropriate. V-band photometry is from the RC3 (de Vau- 
couleurs et al. 1991). For M31, the 2MASS K-band photom- 
etry is likely underestimated (Jarrett et al. 2003) thus the 
V-band total luminosity was used for the stellar mass cal- 
culations. A total K-band estimate for the Milky Way could 
not be found, so the V-band mass estimate was adopted. 
The V-band photometry for NGC 1399 is from McLaughlin 
(1999) for reasons discussed therein. All galaxy photome- 
try is corrected for Galactic extinction. Host baryon mass, 
^baryon, is the sum of M, aud a host gas mass estimate 
that depends on the host V-band magnitude (derived by 
McLaughlin 1999; see Eq. 1 in i)4.4|l . Halo masses, Mhaio, are 
derived using weak-lensing Mhaio estimates (Mandelbaum et 
al. 2006; see E3J. 

The GC total number for NGC 1399 was derived using 
the surface density profile from Bassino et al. (2006B). As- 
suming the GCLF form found for NGC 1399 by Dirsch et al. 
(2003), the Bassino et al. (2006B) surface densities were in- 
creased by 60% to account for incompleteness in their study. 
The best-fit de Vaucouleurs profile was integrated to the GC 
system edge (45'; Bassino et al. 2006B), which yielded a to- 
tal GC population estimate of 5800 ± 700. 

^ http:/ /ncdwww. ipac.caltech.edu/level5/NEDlD/intro. html 



Table A3. Mass-to-light ratios used to compute galaxy masses. 
See text for derivation. 
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V-band 


1.45 
10.23 


1.14 
5.69 


1.10 
5.93 


1.05 
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1.03 
5.23 


0.99 
5.09 


0.96 
5.05 
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0.96 
4.98 
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5900 ± 721'' 


40'' 


N1407 


E 


32. OO'^ 


-25.32 


-22.26 


11.60 
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20.89 ± 0.72 
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1775 ± 1509 


309 


N1399 


E 


31.34" 


-25.03 
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12.25 ± 0.70 
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270 ± 68^ 


30J 


N1052 


E 


31.28" 


-23.85 


-20.84 
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80 ± 20P 


36?' 


N2683 


Sb 


29.44" 


-23.14 


-20.47 


10.58 


10.59 


11.88 


2.48 ± 0.04 
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Table A2. Globular cluster specific frequencies. See Table[7]for parameter definitions. 
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6.25 ±0.76 


3.75 ±0.46 
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8.05 ± 1.15 


4.83 ±0.69 
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3.13 ±0.45 
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M31 


2.77 ±0.62 
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0, 
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0, 
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9.42 ± 1.14 


4.71 ±0.57 
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0, 
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17.73 ±2.15 


8.86 ±1.08 


8.86 ±1.08 
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1, 
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N821 


2.94 ±0.41 
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0, 


.88 ±0.12 


2.79 ±0.39 


1.95 ±0.27 


0, 
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0, 
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1.76 ±0.45 


0, 


.75 ±0.19 


2.40 ±0.61 


1.68 ±0.43 


0, 
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0, 


18 


8.35 ±2.12 


5.84 ±1.49 


2.50 ±0.64 
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3.89 ± 1.17 


1.94 ±0.58 


1, 
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3.73 ± 1.12 
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1, 
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0, 
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13.06 ± 3.92 
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0, 
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1.70 ±0.19 
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0, 
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0. 
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0, 


.91 ±0.15 


4.83 ±0.82 


3.96 ±0.67 


0, 
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0, 


15 


17.32 ± 2.95 


14.20 ±2.42 


3.12 ±0.53 
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4.36 ±0.30 


1.05 ±0.07 


3, 


.31 ±0.23 


4.21 ±0.29 


1.01 ±0.07 


3, 


.20 ± 


0, 
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15.18 ± 1.05 


3.64 ±0.25 


11.53 ± 0.80 
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4.32 ± 0.44 


1.94 ±0.20 


2, 


.37 ±0.24 
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0, 
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18.48 ± 1.89 
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0, 
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0. 
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0, 


.63 ±0.16 
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3.14 ± 1.05 


2.10 ±0.70 
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1, 
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0. 
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